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We present detailed calculations of the moments of structure functions, performed 
explicitly in a covariant gauge starting from familiar techniques. The results illustrate the 
central features of perturbative QCD (universal factorization of mass singularities, lad- 
derization, etc.) and provide insight into their origin. 


1. Introduction 


In the last two years QCD has been shown to support a parton-like picture of 
certain hard scattering processes. The damping of large momentum transfers, due to 
the decrease of the running coupling constant g, is just sufficient to ensure that the 
parton-model diagrams dominate to leading order in g*. Furthermore corrections 
can be calculated in a systematic expansion in powers of 2”. 

Given a large momentum scale Q we can always use the coupling constant Z(Q), 
but this is not enough to justify the use of perturbation theory. Technically a mea- 
surable quantity W generally depends on masses (7), as well as any other dynamical 
variables (Q’) involved, i.e., W= W(g(Q), O/m, O/Q'). If we attempt to expand W in 
powers of g for large Q we may find powers of In(Q/m) which spoil the expansion 
(it may also be spoiled by powers of In(Q/Q’) if Q/Q’ is not held fixed when Q is 
made large). Physically, sensitivity to m indicates that the large-distance behaviour 
of the theory is involved, which we cannot possibly treat perturbatively. Therefore 
the use of perturbation theory requires either that In(Q/m) terms do not occur (an 
example of such a mass finite case is of9') or that the logs can be summed up and 
the sensitivity to m can somehow be factored out. 

This “factorization” occurs in deep inelastic scattering where the operator prod- 
uct expansion (OPE) implies that the moments of the structure functions ** W, 


* Address after September 1: Physics Department, Rockefeller University, New York 10021, 
USA. 
** Tn our explicit calculations we consider W,, = { x"F4(x, Q*) dx =f (x"—-1F,/2) dx to leading 
order in £. 
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satisfy 


N 
OQ? Q? 7 Q? m? m2 
Wale, 25% =2i CP 8,3 M g"9 * Oz, : (1.1) 
where yp is the arbitrary scale chosen to define the coupling constant g. Because y is 
arbitrary * 


dw ) og a 
GE (aya aytag) Man 
In the non-singlet (NS) case (in which N = 1 in (1.1)) this leads immediately to the 
renormalization group equations ** (RGE) for C” and M”, the anomalous dimen- 
sions being separation constants of the differential equation. The RGE allow us to 
sum up the logs of Q?/u? and obtain the behaviour of the C” for large Q? in asymp- 
totically free theories. 

The separation into a calculable Q? dependent factor (C”), determined by the 
short-distance properties of the theory, and an incalculable long-distance factor 
(M") can be put into closer correspondence with the parton picture by introducing 
another large scale M and setting u? = Q?: 


Wn°(Q7)= Chs(8, O7/u?)[CXs(g, M?/u?)1—! WS (M?) 
= CRs(8(Q7)){ [Xs @(Q?), M?/0?)]—* Wr (M?)} . (1.2) 


The first factor (C”) represents the “hard scattering cross section”, while the 
second (C")—! W,, is the “Q? dependent quark distribution”, with calculable Q? 
dependence but unknown magnitude. (In the singlet (V = 2) case we need to con- 
sider simultaneously the function W,,, which would be measured by currents cou- 
pled to gluons; we then obtain coupled RGE and a description in terms of quark and 
gluon distributions.) 

Recent progress has consisted of two closely linked parts: 

(i) analysis of the diagrams summed up by the RGE, leading to information 
ab out the topology of the final states Gets, etc.) which is beyond the scope of the 
OPE; 

(ii) the demonstration that factorization occurs in other hard processes (e.g., 
heavy uw pair production and large-py; hadronic processes), the same “universal’’ 
M" occurring in each case. Likewise in production processes (e*e~ > 7 + ..., etc.) 
there is factorization and the long-distance part is universal. 

In this paper we present calculations in the leading log approximation (LLA) 


* If the quark mass parameters m are defined at uw there should also be a u2(am/apu2) a/a Mg 
term but this contributes to terms of order m2/Q? in (1.1), which we ignore. 
** The relevant equations are collected in the appendix where we also define all the symbols 
employed in the text and figures. 


A.B. Carter, CH. Llewellyn Smith / Perturbative QCD 399 


which support the idea of universal factorization, and, we feel, give insight into the 
structure of the theory. An outline of our approach was published some time ago 
[1]. It turned out that very similar conclusions had been reached independently by 
Dokshiter et al. [2], who (like us) were inspired by the classic work of Gribov and 
Lipatov [3], and by others [4], many of whom were inspired by the important 
papers of Sterman and Weinberg [5], on mass finite quantities, and of Politzer [6], 
who stressed the notion of factorization of mass singularities. Subsequently factori- 
zation has been proved for non-leading logs [7]. Nevertheless we feel that we have 
something to add for several reasons. 

(i) We work in a covariant gauge in which we can calculate explicitly. In contrast, 
other authors have mostly employed general arguments and used axial gauges, which 
makes it hard to check the conclusions. In fact, ours is the first fully explicit calcu- 
lation of W,, beyond O(g?) in any gauge as far as we know. (The pioneering work of 
Gribov and Lipatov was done in the Feynman gauge, but they considered an abelian 
model which is much less complicated than QCD; furthermore we believe that we 
are not alone in having found their paper very hard to understand.) 

(ii) Our calculations start from familiar diagrammatic techniques and therefore 
offer considerable heuristic value. The physics is made clear: one sees the infrared 
and “‘Kinoshita-Lee-Nauenberg”’ (KLN) cancellations [8] occur line by line and the 
workings of gauge invariance to eliminate singularities associated with spurious 
(gauge) degrees of freedom. 

We consider as a paradigm the specific case of W,, for deep inelastic scattering (it 
should be clear that the technique can easily be applied to other hard processes). 
Our goal is to rederive the well-known result that 


= NS 492 
m0) -[Rom)  MSOr) (13) 


in LLA and the corresponding singlet result (the symbols are all defined in the 
appendix). For a given power (g?)" we identify the leading powers of In(Q?/u?), 
In(Q?/m?) and In(m?/y?) and then sum over n (some terms can be constructed 

from others, using the fact that W, must be y-independent, but we prefer to identify 
their origin explicitly). 

There are four possible sources of logs: 

(i) ultraviolet logs, due to the dressing of vertices and propagators and to wave- 
function renormalization (these logs introduce u-dependence); 

(ii) logs due to collinear radiation of real particles (which depend on m); 

(iii) ‘parallel logs’ due to virtual collinear emission and absorption of almost real 
particles (these also depend on m); 

(iv) logs left over after the cancellation of infrared divergences between real and 
virtual processes involving soft particles (in fact we shall show that such logs are 
absent). 

If we consider the virtual forward Compton amplitude 7 (recall that structure 


400 A.B. Carter, CH. Llewellyn Smith | Perturbative QCD 


functions W « Im 7), it turns out that intermediate N-particle states in the t-channel 
only give rise to mass singularities (giving logs of types (ii) and (iii)) if NW — 2 of the 
particles are longitudinal gluons, with €,,(k) «k, (this is easily made plausible by 
dimensional arguments [7]; we shall see it explicitly). It is clear that the contribu- 
tions of these longitudinal gluons are in some sense artificial. Mass singularities are 
associated with the propagation of real intermediate states [9], i.e., virtual processes 
which could become real as m > 0, allowing the vertices to separate by arbitrary dis- 
tances. Longitudinal gluons are unphysical and cannot contribute to genuine real 
processes. Therefore, if we sum over all diagrams, their contributions must cancel. 

In the axial gauge the propagators of longitudinal gluons have no poles as k* > 0 
(and correspondingly there are no ghosts). Therefore, only states with V = 2 give 
singularities diagram by diagram (‘“‘ladderization’’). The price to be paid for this sim- 
plification is that there are additional (spurious) poles in the propagators and explicit 
calculations are very difficult. 

Insofar as real radiation is concerned, we can reproduce this simplification (regard- 
less of our choice of gauge) merely by summing over physical, transverse polarization 
states for real gluons. We sum up the effects of virtual (longitudinal) gluons using 
Ward identities, which express the fact that contributions of longitudinal gluons 
must be gauge artifacts, since they are not associated with a physical pole. 

Covariant gauges are generally characterized by a parameter a such that the gluon 
propagator Dy, satisfies 

k? Dw(k?)| k2s—p2 > CS ww t (= 0) (Su 7 a 
It turns out that (as discussed in detail in sect. 3) the k,k,/k? term gives rise to 
problems in our approach. They cannot be avoided by choosing a = | since the kyk, 
term is induced by radiative corrections for k? # —u?. However, it is possible [10] 
to let a be a function of k? such that 


k? Diy (Kk?) = —0(k?) Sup - (1.4) 


We call this choice “stagnant Feynman gauge” (SFG). In SFG the usual renormaliza- 
tion group equations hold, but without the “gauge drift” term 60/da and with 6(g) 
and the y(g)’s calculated in SFG (this is the same as calculating them in the Feynman 
gauge to the order considered in this paper). 

This paper is organised as follows. In sect. 2 we discuss the calculation of W,, to 
O(g’) in detail in order to develop calculational techniques which can be used in 
higher orders. In sect. 3 we calculate WNS to O(g*); this is enough to understand 
essentially all the physics of asymptotic freedom and mass singularities. In sect. 4 
we present the Ward identities and use them to sum up the effects of “parallel logs’ 
to all orders. The summation of all contributions to all orders (in LLA) and factori- 
zation using a large scale M are discussed in sect. 5. Sect. 6 is devoted to mixing and 
the singlet case. In sect. 7 we comment on the extension to non-leading logs and 
applications to other process. 
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2. Order g 


In this section we discuss the leading logarithmic contributions to W,, in O(g?) 
introducing techniques needed in higher orders. Although the non-abelian couplings 
are not yet evident, much of the essential physics appears already in this order. At 
the end we shall display the features which persist in higher orders. Initially we shall 
encounter infrared singularities as formally divergent integrals. We discuss the regu- 
larization of these integrals and their combination into an infrared-finite result when 
all the ingredients are assembled. 

With m = 0 the amplitudes of figs. 1a,b would be singular as tf > 0 and ass >0 
(x = Q?/(2p - q) > 1) respectively, since at these points the exchanged quarks 
become real and the propagators blow up. Provided we only consider physical 
gluons with transverse polarizations, the vertex in fig. 1a is proportional to /—f, 
because the quark-gluon vertex conserves helicity (to order m) and emission of a 
transverse gluon is forbidden as 86 > 0 (t > Q). This means that for x #1 the diver- 
gence in W comes from the square of fig. 1a only and is logarithmic (the additional 
singularity as x > 1 coincides with an infrared divergence as discussed below). We 
can find the coefficient of this singularity by extracting the (—t)!/? factor from the 
amplitude and ignoring all terms of higher power in ¢. Thus, dropping all terms 
which do not contribute to the logarithmic singularity as m —> 0, the squared matrix 
element (summed/averaged over spins and colours) is 
(P — F) 40 PB ¢E 


My=(ig)??Ce 22 Trp 
a (ig) Ehy eae  deeair errr —ttm? 


= (6)? Ce ooa) D Tele vy @ ~ B= Hb + 0-97) 


2 1+ty? 
~ 2C 
BNF “+m? yl —y) 


Tr(p — §) M°, (2.1) 


where (i) we assume Q* >> m? and work in the frame q° = 0 in which the momenta 
are large: 2E ~ 2|p| = /Q/x; (ii) we can therefore set ky = (1 —y) Py + O(/2) for 
y #1; (iii) the gluon polarizations are spatial and transverse: e? =0,€-p= 

[pl sin @ ~ —t/|Kl; Civ) Tr(p — &) Mo is the squared matrix element for a real incom- 


Fig. 1. 
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ing quark with momentum fraction 
’ Q° ed 
2q-(p—k) y 


in our leading log approximation. In this configuration the gluon phase space simpli- 
fies to 


dk 
2k Qn = 


5 dyd(-t), 


so that finally we obtain 


pr, ity 4. 
W(x) =ACpln = f dy——=— w(x’), (2.2) 
Spe ed say) 


with A = g?/8n, where the log arises from the ¢ integral. Taking moments, 


mes fs x"W,(x) dx =) Cp In => fort “ely we. (2.3) 


The Perr as y > I reflects the infrared divergence as |k| > 0 and will be dealt 
with below. Eq. (2.2) clearly generalises to the more complicated process shown in 
fig. 2, with W(x’) replaced by the probability corresponding to the process quark + 
photon > X induced by a quark with momentum yp. . 

Similarly, we can extract the log contributed by the square of the amplitude of 
fig. 1b by noting that in our frame, as s > m? (x > 1), p’ becomes parallel to k, so 
that we can set kK, =(1 —y)(@y + Pu) after we have extracted a V/s “helicity conser- 
vation” factor. Thus fig. 1b ate 


ee ny + BE bey. (2.4) 


My x g? Cre — we t= 


Introducing the identity 
1=84(pt+q—) dl &(2 —s) ds, 
the phase-space factor is proportional to 


54(p+q—1) d4l (2 — s) dvds 


Fig. 2. 
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Pp’ 
k 
X 
P 
Fig. 3. 
in the p’ / k configuration. Hence we obtain 
ds 1lty?f, (s — m*)x 
Wy (x) =A Ce f dy [——4 hey [pou - x8 : (2.5) 


where the last factor would be just Wo(x) for s = m?. In fact, since we only want 
the x-moments, we can replace (2.5) by 


1+ y? 
1 — 


4 Wo(x) (2.6) 


OT fp 
Moo) =A Orin Fa fw 


in LLA if we realize that the s-integral converges for y # 1 since ys > (1 — y) m? 
(the infrared divergence is discussed below). It should be clear that in LLA the con- 
tribution of the process of fig. 3 is given by eq. (2.6) with Wo(x) replaced by the 
probability corresponding to the process quark + photon > X + quark and Q? 
replaced in the logarithm by some suitable scale set by the process within the circle. 

The interference of amplitudes 1a and 1b does not contribute to leading logs if 
the emitted gluon is transverse *. The implications of this important point are dis- 
cussed further below. 

Fig. 4 contains both mass singularities and an ultraviolet divergence. Since we 
perform the vertex renormalization off-mass-shell, the two kinds of divergence are 
easy to separate. The mass singularities correspond to real subprocesses, where inter- 
nal particles propagate on their mass-shells and the vertices can separate to arbitrari- 
ly large distances in coordinate space in the m = 0 limit. For fig. 4 there are two dis- 
tinct possibilities (modulo the infrared region which we discuss below): k? = 0 with 
k® > 0 and either (i) (p — k)? ~ 0, p® — k° > 0 or (ii) (p’ — k)* ~ 0, p’® — k9 > 0. 
In order to see this, parametrize the loop as follows **: 


Ky =z(xPputaqy) t+ (1 —y)Patki, (2.7) 


* If we included the (spurious) contributions of scalar and longitudinal gluons (which cancel) 
by using Le wey = ~fpp> the helicity argument would break down and the interference term 
between la and 1b would also give a log. 

** The choice of the projection vector given by the first term on the right-hand side of (2.7) is 
more or less arbitrary. One merely requires that its square be of O(m?) and its contraction 
with py be large and positive. 
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harks 


Fig. 4. 


choosing kt, such that kt - p =k'- (xp +q) =O and k'- kt =—k?. Then focusing upon 
the propagators contributing to possibility (i) above and setting m = 0 for the mo- 
ment, 


d*k 7 2nv dzd(1 — y) dk? 
(k? +ieX(p — k)? tie) (2vz(1 — y) — k? + ie)(Qvzy — k} + fe) ’ 


where v = p -q. The integral is zero unless the integrand has poles in both the upper 
half complex z-plane and the lower half z-plane; thus y and y — 1 have opposite 
signs, i.e., 0 <y <1, corresponding to k® > 0, p® — k° >0. No matter which plane 
we close the integration contour in, corresponding to putting one or the other propa- 
gator on-shell, the residue is proportional to 1/k?, and the log comes from integra- 
tion of this variable over an asymptotic range bounded above by Q? and below (rein- 
stating the quark mass) by m?. 

Since z ~ k? in the vicinity of either pole, all terms of O(z) or O(K?). are neglect- 
able in LLA; in the numerator, we can set k,, =(1 — y) py. Thus at the incoming end 
of the numerator spinor chain we can write 


(2.8) 


uo G—f- m)yyu(p) = ok uu(p). (2.9) 


This means that the virtual gluons responsible for the logarithmic divergence are 
longitudinally polarized (€, ~ k,). When this gluon is attached to the other end of 
the spinor chain, we get 


' 1 , 1 
WP) Yp a KE (p(B — m) — @ - F - m)] 
“p-—k-m p-k-m 
=-a(p')~, 
so that to LLA the amplitude is proportional to the Born amplitude. This is the 


first primitive example of the use of a Ward identity to treat longitudinal gluons. 
Finally, then, the contribution of fig. 4 corresponding to k / p (the first possibili- 


ty) is: 


wir f aa for | w= rcein Sl f ay] we, 
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An identical contribution is similarly obtained from the k / p’ region: 
KI a 
WkIP' = Ce n& fo = we. (2.10) 


We pause now to consider the problem of soft (infrared) divergences. Adding 
eqs. (2.3), (2.9) and (2.6), (2.10), we see that the divergences as y > 1 cancel sepa- 
rately in the k / p and kK / p’ regions, and this naive procedure actually gives the cor- 
rect answer: 


2 
W23*2-9) =) Coin 2 ele ae ie 02) — 2) yp 


= Cr(f, + 5) In Zw, (2.11) 


1 


: 1 + y?) — 2y 
Wh2-6)* (2.10) =rC Q d area 0 
F In m2 f 4 1 =) Wr 


1, @? 
=i Cr) In we, (2.12) 


(where f,,, which is defined in the appendix, is proportional to the familiar non- 
singlet anomalous dimension of the operator product expansion). We shall now show 
that when we combine real and virtual contributions the region k < 5/Q? (k =|k\) 
contributes terms of order 5 In(Q?/m?) (where 4 is arbitrarily small but fixed as 

Q* >), We can therefore exclude the region y > 1 — 6 and the region where k is 
small but @ is large (which we have ignored altogether). It follows that (2.11) and 
(2.12) are correct. 

Infrared divergences are most easily treated by (temporarily) using Le,€, = —S yp, 
in which case it is well-known that there is a direct topological correspondence 
between real and virtual contributions which combine to give finite answers. Treat- 
ing them in turn: 

(i) | fig. 1b |? + (p’ wave-function renormalization to elastic scattering) is finite. 
The divergent contribution of fig. 1b is 

Djs) fot 

ee’ k \p'-k k - (E'—p' coséy ~ 
The first factor behaves as (Q7)‘/e if we go to 4 + 2e dimensions (€ > 0). The angu- 
lar integral would be linearly divergent if 7m were zero, so the second factor is pro- 
portional to m?(1/m?)(m?/Q?)©. The net result (~(m?)£/e) depends only on m?, as 
we might expect for incoherent radiation from an external line, and cancels the 
divergent contribution to wave-function renormalization. 


406 A.B. Carter, C.H. Llewellyn Smith / Perturbative QCD 


(ii) Likewise | fig. 1a|? + (p wave-function renormalization to elastic scattering) 
is finite and does not leave a In(Q?/m7?) factor as might have been feared. 

(iii) The interference of figs. 1a and 1b + fig. 9 is finite. As usual the k and angu- 
lar integrals factorize, the former giving (Q7)*/e. The angular integral is at most loga- 
rithmically divergent as m — 0 so that in both the real (r) and virtual (v) cases it 


gives 
Br’ m € 
AR 7) - 1 ; 213 
alt ay 
but it cannot generate terms ~(m?/Q7?)€ (which could only come from m? times an 
integral which is linearly divergent as m > 0). The sum of real and virtual contribu- 
tions is infrared finite, which requires that separately 


BY + BY=0 
(to eliminate terms ~e—! In(Q?/m7?)) and 
At +AY=0 


(to eliminate e—!), with the result that no In(Q?/m?) terms remain. 

It follows that the region k < &/Q? contributes terms of O(5) + O(6 In(Q2/m?)) 
(the latter being due to the usual parallel divergences unaccompanied by k—!) and 
can be ignored. The same is true in the non-abelian case, including the effects of 
radiation from gluon lines, because real and virtual contributions still cancel between 
diagrams with related topologies [11]. Again wave-function renormalization plus 
incoherent radiation from external lines is finite and the cancellation does not leave 
logs since both behave as (mm?) e—!. Virtual contributions combine with related 
interference terms to a finite result free of logs because both behave as (Q2)€ e—! 
for the reason discussed above. We can therefore freely manipulate divergent y-inte- 
grals, which will combine into an infrared-finite result, and ignore the k small, é 
large region, as we have done above (thus it suffices to consider all gluons to be 
“hard’’), 

Finally, we must include in W,, the “ultraviolet” logarithm from the k > © region 
of the loop in fig. 4, which gives 


2 
WUY =) Cr ln (5) we (2.14) 


and the ultraviolet divergent contribution * of wave-function renormalization on 
each of the outgoing lines, given by the usual renormalization group anomalous 
dimension: 
2 
wren =9 Er Ge in“) | we. (2.15) 


* We ignore the infrared divergence in the renormalizations. They cancel against soft diver- 
gences in |1a|? and |1b|”, as discussed above. 
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At last we add up all the contributions (eqs. (2.11), (2.12), (2.14), (2.15)) to W,, in 
this order: 
2 


WOE) = Cp fy in(S) We =n, in( 2) Wr (2.16) 


To O(g?) accuracy, then 


(nsf) rasa) atsouy, en 


which is written in the (trivially) factorized form of eq. (1.1) and satisfies to this 
order the appropriate renormalization group equation. 

We now call attention to several important features of this calculation which per- 
sist in higher orders. 

(i) With only physical (transverse) produced gluons interference terms such as 
those illustrated in fig. 5 never produce leading logs. This follows immediately from 
the “helicity” zeros which make the individual amplitudes behave as Vs/s, Jit etc., 
when we approach the parallel configuration. Since these factors diverge in different 
regions of phase space, only symmetric diagrams (as in fig. 6) possess the required 
logarithmic singularities. The fact that interference terms do not contribute in LLA 
(“‘ladderization’’) allows a probability interpretation and accounts for the success of 
the parton model. 

(ii) The leading logarithm comes from a wide (asymptotic) range of the logarith- 
mic variable (s, t, etc.). We can truncate the integrals at either the upper or lower 
end and still obtain the leading log asymptotically, e.g., 


eQ? 


2 aittee © altlut 2 alt! 
LL 
f ee = i —————_- = i, vitae ts (2.18) 
Jo Ith Jo Iel go, Te 


where € is arbitrarily small but fixed as Q? > c, Nevertheless, the transverse momen- 
tum (p?) ~ <|t|) vanishes to this order of approximation, since the integral (2.18) is 
no longer logarithmic if any power of ¢ is added to the integrand. In this sense, it is 
the region of small logarithmic variables which gives the leading logs. This means all 
particles are emitted at small angles to the vectors p and p + q (@ > 0 in fig. 1), and 


et kf Rf 
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Fig. 6. 


the final state has a two-jet structure. Thus in LLA we can treat the particles 
involved as highly virtual (|¢| >> m7), and use the asymptotic form for ultraviolet 
vertex and propagator corrections, yet almost real (|t| << Q*) as far as kinematics 
is concerned. 

(iii) The virtual gluons contributing parallel divergences are longitudinally polar- 
ized (€, ~ ky), and the amplitudes contributing these singularities can be summed 
up using Ward identities. The result is an “effective” emission diagram (see fig. 7) in 
which the divergence is collected on the external line. With this effective diagram is 
associated a factor 


1 
r ce( f ay) X log , (2.19) 
0 


as in (2.9) and (2.10). 
(iv) Mass singularities associated with the outgoing quark (eqs. (2.6), (2.10), 
(2.15)) cancel in the inclusive sum, leaving a residual logarithm 


Q? 
5 Cr in( 2) (2.20) 
a 
This cancellation is the realization, in this order, of the Kinoshita-Lee-Nauenberg 
(KLN) theorem. The log, (2.20), which remains is the ultraviolet renormalization 
appropriate to a quark propagating an amount Q? from its mass-shell. This is a gen- 
eral feature: there are never any mass divergences associated with outgoing (i.e., 
unobserved) lines, and the In(:?/m*) terms associated with wave-function renormali- 


Ward 
> 
identity 
insertions Ey ku 


Fig. 7. 
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zation get turned into In(u?/M7), where M is the upper range set by the kinematics 
for the mass of the line and all particles emitted from it (in the LLA, there is no 
other scale). It is plausible (and we shall later show in detail) that the total effect of 
radiation (real and virtual) off unobserved lines, is given by the functions Z3 and Z§ 
which multiply the bare propagators for very virtual massless quarks and gluons. 
(The functions Z3 and Z£ can be calculated using the renormalization group). 

(v) The soft region (k < 5/07) can be ignored completely. We can therefore 
freely manipulate apparently divergent y integrals which combine to well-defined 
results, e.g., the k / p logs give 


1 + x? * [l+y? 1 
58-39 fa] -5. 


1 — l—-y 2 


which is the Altarelli-Parisi distribution [12] apart from the final factor (—4). In 
higher orders it is always the finite moments of this combination which occur in 
W,,, whose n dependence is determined from f, ~ YN), as in eq. (2.11). 


3. Order g* 


In this section we present explicitly the leading logs contributed by all O(g*) dia- 
grams. This is enough to understand essentially all the physics of asymptotic free- 
dom and mass singularities in QCD perturbation theory. After some preliminary 
remarks we shall] discuss all of the various contributions (ultraviolet or parallel) of 
every diagram in turn. The results are given in a table enabling us to focus on the 
physics in the text. At various junctures, which we denote CC (for “collection and 
comment’’), we shall pause to collect various contributions which combine together 
in a significant way, with an eye to features which generalize in higher orders. The 
complete result is assembled and discussed at the end. 

Up to now we have treated massive on-shell quarks, but henceforth we set m = 0, 
taking incident quarks (and gluons, in the mixing case) slightly off-shell. It is also 
convenient to take the final particles slightly off-shell *. This is largely a bookkeeping 
device, regularizing and labelling parallel] divergences at intermediate stages of the 
calculation and indicating how the KLN cancellations work. Since in reality (as dis- 
cussed in sect. 5) we are provided with virtual initial quarks and gluons by a hadronic 
vertex which we cannot calculate, and in any case the notion of mass-shell is not sen- 
sible for confined quarks and gluons, it is really the structure functions for virtual 
particles which we require. Of course, off-shell quantities are not well-defined (e.g., 
they are gauge dependent). One is accustomed to thinking that the difference 


* In principle the incident particles should be spacelike and the produced particles timelike to 
avoid spurious singularities. In practice these cause no trouble and we shall not bother to 
keep track of signs and distinguish, e.g., In(Q2/p2) from In(Q?/—p?). 
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Fig. 8. 


between on- and off-shell quantities must be of order (p? — m*)/Q? and is therefore 
irrelevant at large Q?. Happily, in our stagnant Feynman gauge this is true diagram 
by diagram, a circumstance which would not obtain in a general covariant gauge *. 
We can therefore treat diagrams with off-shell particles, confident that our calcula- 
tions continue smoothly onto the mass-shell, where their gauge invariance is assured. 

Turning now to explicit evaluation of the O(g*) diagrams, we begin with the 
instructive and relatively simple case of qq production (fig. 8). These diagrams give 
double logs in the configuration p/p ((p + P')* =S ~~ 0) and p +p’ parallel to p or 
p' (t ~ 0 ors ~ 0). The amplitudes behave as s~!/? times ¢—!/? or s~1/? (from heli- 
city conservation) so that (as in the case of figs. 1a and b), 

CCI: we need not consider interference terms to obtain leading logs. Because the 
qq vertex conserves helicity the virtual gluon becomes transverse as p and p’ become 
parallel. Thus it is clear that the contributions of figs. 8a and b are proportional to 
la and b, respectively. 

CC2: In diagrams with the topology of fig. 8 (e.g., fig. 11 below) the singularities 
as § > O are associated with the possibility of virtual processes becoming real. The 


* In other covariant gauges the Kyky|(k2)? term in the propagator creates problems. These are 
illustrated by the amplitude of fig. 4. If we set = (f — m) — (pg — & — m) at the bottom 
vertex, the first term is equal to zero (using the Dirac equation) times an integral which is not 
only linearly divergent on-mass-shell but also contains a parallel divergence, iLe., we find 
(0/0)(@ + b in(Q2/m2))! It is standard practice (e.g., in proving the equivalence of different 
gauges) to drop this term, ie., it is implicitly assumed that the theory is provided with an 
infrared cut-off which makes such ambiguous terms equal to zero on-mass-shell (e.g., by 
going to 4 + € dimensions, in which case 0/0 > (p? — m2)!t€/(p2 — m2), or by giving the 
gluon a mass A, in which case 0/0 > (p2 — m*)/(p2 —m? + 42)). However, if we regulate the 
infrared divergences simply by going off-mass-shell, we will clearly obtain the “wrong” answer 
when we return to the mass-shell (0/0 > (p2 - m?)/(p? — m?) = 1). Presumably we would get the 
right answer for infrared-finite quantities (e.g., o(y"q >q)+ o(y"q —> q + soft gluon)), in 
which the translation of # to produce # — m can introduce at most a logarithmic divergence, 
provided we treat real and virtual particles in the same way. Thus in gauges with kyk,/(k)? 
terms in the propagator we must either abandon our decision to consider only the produc- 
tion of transverse gluons and learn how to handle the discontinuity of kykyl(k?)? or stay 
on-shell (alternatively we may finesse the problem by interpreting (p2 — m?)/(p2 — m?) as 
zero off-shell when we anticipate that such terms would vanish or cancel on -shell {13]). 
These problems occur also in the Feynman gauge in O(*) and beyond when Kuk, terms are 
induced in the propagator by radiative corrections. They are avoided altogether in the stag- 
nant Feynman gauge, in which these terms are absent diagram by diagram. 


A.B. Carter, CH. Llewellyn Smith / Perturbative QCD 411 


Fig. 9. 


coefficient of the singularity is therefore proportional to figs. 1a,b evaluated with 
real physically polarized gluons. 

Doing the (p — p ) phase space integral we find that the constant of proportion- 
ality is 


2) Tr In(M?/p?) , 


where p* = p’? is the mass of the produced quarks and M? is the invariant which 
controls the upper range of Ss: ¢ in the case of 8a and s in 8b. 

CC3: These contributions should be combined with the qq contribution to the 
wave-function renormalization of the produced gluon in fig. 1, indicated symboli- 
cally in fig. 9. Renormalization produces aN Tr In(p*/y), so the sum contains 
In(M?/?) only and is ‘“‘mass finite” as 5* > 0 (KLN theorem). The net effect of qq 
radiation from the gluon line of fig. 1 is thus to modify it by the qq contribution to 
the renormalization factor Z3(M*/u*) appropriate to a gluon propagating a distance 
M? from its mass-shell. We shall see that radiation from outgoing lines always con- 
spires to produce such an “effective” off-shell propagation factor. 

Integrating over ¢ (or s), we obtain (we shall quote all results in units of A? W®): 


2 


2 2] } 
gat9a — 1 22° 2F| nity 
ws kCeTe [In iar a are (3.1) 


with the same contribution from 8b and 9b with p? > p’? and n = 0 in the y integral. 
Next consider the two-gluon production diagrams of fig. 10. Simple kinematics 
establishes that in 10a and 10c: 
CC4: —t' >—t,s >s’, so that the integrals over these variables are nested, e.g., 


Fig. 10. 
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This is the beginning of an exponential series (1/n!)(log)” for such diagrams. 

The resulting contributions of 10a and 10c are given explicitly in table 1. 

Two gluons may also be produced through the gauge coupling in this order. This 
process (fig. 11) gives double logs when k / k’ ((k +k’)? =3 ~ 0) andk+k' /porp' 
(t ~ 0 or s ~ 0). The individual amplitudes behave as §~1/? times ¢—1/? or s~'/?, so 
that each diagram contributes incoherently to leading logs (CC1). The yy at the ver- 
tex indicated must couple to €,, Eu Ky, Or Ky. Consider these factors in the double 
log configuration. Both k,, and k,, are proportional to p, (or Dy) and may therefore 
be dropped using the Dirac equation. Since ¢ and &’ are perpendicular to k + k’, it is 
clear that the contributions will be proportional to those of fig. 1 (CC2); integrating 
out the k — k’ phase space we find that the constant of proportionality is 


C nt fo Feexeeees (i | (3.3) 

G k?2 J y Wie y y y ¥ ’ . 
where k? = k’? is the mass of the produced gluons, k,, ~ (y/(1 — y)) Ky in the paral- 
lel configuration, and (as in the case of fig. 8) M? is t or s. Shortly we shall find the 
contribution which renders the y-integral finite (CCS). Then (CC6) we shall see that 
the k? dependence is cancelled when we consider further O(g) renormalizations of 
the gluon wave function. 

Next we consider virtual radiative corrections to fig. 1 shown in fig. 12. First 
there are the ultraviolet (|k| > %) logs due to the loops in all but 12e,j,k. These 
simply modify the O(g?) emission amplitudes by factors of In(t/u?) (12a—c), 

In(s/u?) (12f~-h) or In(Q?/u?) (12d, i). Integrating over s or ¢ they give the contribu- 
tions listed in the table. 

These ultraviolet modifications are summed up in the shaded circles in figs. 13a—f; 
in this figure and henceforth in similar diagrams, shaded circles stand for ultraviolet 
corrections only. 

The diagrams of fig. 12 also contribute to leading logs when the virtual gluons 
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Fig. 12. 


are near their mass-shells and emitted parallel to the quarks (real gluons are under- 
stood to be parallel to the line from which they are emitted). Let us consider dia- 
grams 12a,b,d,e,i,j,k in the configuration k2 ~ 0, k / p. (By the p * p’ symmetry 
of fig. 12 it should be clear from this discussion how to compute the contribution 
of radiation parallel to p’.) In this configuration the virtual loop in 12d gives a 
In(Q?2/t), to which we associate the effective diagram 13g (see fig. 7 for the notation 
used in fig. 13). We can evaluate this log and do the integral over ¢ immediately 
using (2.2) and (2.19) to obtain 


W138 = 3 CE In?(Q?2/p? rn Kv) , (3.4) 


DET ET 
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where 
1 
—2y 
v=] dy-—. 


As in the case of fig. 4, the virtual logs in 12a,b,e,i,j,k arise when k* and (p — k)? 

are small and, since k / p, we may write at the gluon vertex 
2y 

“@-#) yup) = fay tt . (3.5) 
The contributions of these diagrams can then be combined using a Ward identity. 
We shall study these Ward identity manipulations in a more general way in sect. 4; 
in this order our results can easily be checked explicitly. As in familiar QED, the 
simplest QCD Ward identity states that 0,,A¥ decouples from physical states (on- 
shell and transversely polarized in the case of gluons). Our procedure, applied to 
12a,b,e,i,j,k, inserts a 0 - A (i.e., a gluon with €, > k,,) everywhere in the complete 
O(g?) amplitude for y*q > gg with only the initial quark unphysical. It is therefore 
no surprise that the Ward identity produces the amplitude for y*q > qg apart from 
a factor associated with the initial quark (as in the simpler case of fig. 4). The result 
can be expressed by figs. 13h,i, which the reader should now be able to evaluate. 

Figs. 13j,k,] summarize the effects of virtual radiation parallel to p’. 

We have not quite exhausted the parallel divergences due to virtual radiation in 
fig. 12. Figs. 12b,k,g are also logarithmic in the configuration k*, k'2 ~0,k/k’. 
Dropping terms which vanish in this configuration, the trigluon vertex is propor- 
tional to 

(2k + k')y€g — (2k' + kK)g€q & =) ka €g + (2 


) kge, 


in the notation of fig. 14, where we have written ky ~ (y/(1 —y)) ky. The presence 
of the factors ky and kp means that once again 0 - A is involved and the amplitudes 
combine in a simple way due to Ward identities (we again defer further discussion; 
the result is easily checked explicitly in this order). The factor €,,g ensures that the 
result is proportional to the amplitudes for y*q > qg. We represent it by the dia- 
grams of figs. 15a,b, where the dotted line means that a factor 


Mp, (2-y te) £ 
= — pe 3.6 
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B k 
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with M? = t, in 15a, M* =s in 15b, multiplies the amplitude. 
CC5: Combining figs. 15a,b (eq. (3.6)) and fig. 11 (eq. (3.3)) we get a 'y"q > ag 
term multiplied by the infrared finite factor: 


figs. 11 + 15a,b > —2Cg In(M?/k?) . 


CC6: Combining this with the wave-function renormalizations symbolized in 
figs. 15c,d we get a result “mass finite” as k* > 0: 


figs. 11 + 15> —2Cg In(M?/y?) . 


This factor amounts to the “renormalization” (given by Z3(M?/y)) appropriate to 
a gluon propagating a distance M? from its mass-shell, a feature we have already 
noted (CC3) for qq production. 

CC7: The contributions of figs. 8a, 9a, lla, 15a and 15c may be represented by 
fig. 16a, where the blob means multiply the bare diagram by the ultraviolet correc- 
tion to a gluon propagating with invariant mass ¢t. This is an element which will per- 
sist to all orders. (Figs. 8,9, 11, 15b and 15d may be similarly represented.) 

CC8: When fig. 16a and figs. 13a,b are taken together, the result is the same as 
what one would get by calculating fig. 1a with the running coupling constant Z(t), 
where 


2 
L + 3Bo8? In(t/u) 


This combination is, of course, equal to the total contribution of the UV blobs in 
fig. 16b. 


z(t) ~ g? — 5r84 (4 Co — 37 F] in(t/y?) . (3.7) 


Fig. 16. 
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Finally we come to the O(g*) virtual corrections to elastic quark scattering 
(fig. 17). These diagrams contain double ultraviolet (UV) divergences, mixed UV and 
parallel divergences, and double parallel divergences. The “double UV” logs are easi- 
ly evaluated using the renormalization group with the result shown in the table. 

Next consider mixed parallel and UV logs. The UV logs due to superficially diver- 
gent subdiagrams are calculated first; we represent the consequent logarithmic modi- 
fications to the bare vertices and propagators by the blobs in fig. 18. The remaining 
parallel logs can now be evaluated as in the case of fig. 4 in sect. 2. Consider the 
k / p region (the result for k / p’ is easily obtained with p? > p’?). In 18c¢ the UV 
vertex gives In(Q?/u?) and the result is just to multiply the O(g?) parallel log (eq. 
(2.19)) by the O(g?) UV log (2.14). In 18a,b,d the scale of the UV logs is k?/u? 
(where ky was defined in eq. (2.7)). In 18e,f the scale is proportional to k - p'/p? ~ 
Q?/u?. We shall write In(Q?/u?) = In(k?/y?) + In(Q?/k?): the former log goes along 
with 18a,b and d; the latter “remainder” term will be combined with double paral- 
lel logs. The contribution to W,, is obtained by logarithmic integration over ke and 
gives 


2 2 
wiseoittert(e dpy= 10. cea 27 el (i (S 2 s (25) v, (3.8) 
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Fig. 18. 
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with “‘remainder”’ 
w18° f(remainder) = —}Cp Cg In?(Q?/p?) v. (3.9) 


CC9: Eq. (3.8) is exactly the same as the term obtained by evaluating eq. (2.10) 
with g? replaced by the running coupling constant 2?(k?). In fact, we shall see that 
with g > g, eq. (2.10) gives the correct parallel log factor to all orders. 

Last we consider “double parallel” logs from fig. 17 (17a,b,e and h do not contri- 
bute such terms), which we combine with “remainder” term (3.9) above. As with 
fig. 12, the various terms from fig. 17 can be combined using Ward identities. These 
identities will be discussed in sect. 4; in this order it is possible to verify the result 
explicitly, for which purpose it is convenient to proceed in two stages. First take 
k4 pork /p’ in 17c,4d,f,g,i,j. The result of summing insertions of the longitudinally 
polarized virtual gluons and combining the sum with “remainder” (3.9) is to pro- 
duce terms which we represent in fig. 19 (the meaning of these figures should now 
be clear). The case in which the inner loop gives a UV log has already been discussed 
(fig. 18c). The result of extracting the parallel (both to p and p’) logs from the inner 
loop is represented in figs. 20a—c (20d follows rather obviously from 17k). In figs. 
20a and c the logarithmic integrals are rested (cf, figs. 10a,c), e.g., 


Q2 12 kf? 4,2 2 
dk t dk 1 2 
: i i = ( z) (3.10) 
2 


fig. 10a> f — —= == 
k’2 2 
p? p 


ke 2 

which again points to an exponential series. 

It might appear that 20b emerges from these Ward identity manipulations twice 
— once from fig. 19a and once from fig. 19b — but this is illusory. The confusion is 
connected with a choice we have, in using Ward identity manipulations to isolate 
the In(Q?/p?) in(Q?/p’”) contribution of fig. 17, whether to begin applying the 
identity first to gluons parallel to p or to those parallel to p’. Of course the result 
must not depend on this choice — or our effective emission notation of fig. 20 fails! 
— and indeed it does not, but to see this in detail one must realize that only fig. 17) 
contributes In(p?) In(p’?) (when k / p, k' / p', not when k / p’, k' / p). The single 
diagram 20b represents this contribution. Fig. 20b is thus unambiguous and counts 
leading logs correctly. Verifying this statement (which corresponds to independence 


Fig. 19. 
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of the choice described above) in higher orders is elementary but lengthy, and we 
shall henceforth use notation like 20b freely *. 

CC10: Combined with figs 18a,b,d,e, and f, fig. 20 gives the O(g*) term in the 
expansion of the exponential product 


Q? 42 Q* 412 
dk dk 
exp E Cru f B 7) exp]. Cru . @ 7K) : (3.11) 
p? p'? 


These exponentials will play an important part in the final result; we shall denote 
multiplication by these exponentials by the notation of fig. 21; a hybrid of our 
effective virtual emission notation [indicating v = i dy(—2y/(1 — y))] and our 
ultraviolet blobs (indicating 2”). 

We have now evaluated all the O(g*) contributions. It remains to assemble them 
all into a final result. First let us show that all singularities associated with p’? can- 
cel, as they should from the KLN theorem. The double p’”? logs are collected in 
fig. 22. They give the infrared-finite result 


Wr? = Leetv + To) Ino 5) “2D si CF n( m3 7)) 


Next evaluate fig. 23, where 23a is analogous to fig. 16b, and 23b (from which the 
remainder is dropped) contains 18a,b,d,e,f in the k 7p’ configuration, giving 


Wa? = —4l 3 Cp — 3Cg] [Cro tr0)] in? (S) - m() 


=—1¢, [4 Cus $Ca]| i? (2) - in(25)] 


Last, evaluate figs. 24a and b, corresponding to 13e and 18c (with k / p’) respec- 


* The proof is a generalization of the statement, made above, that 17j does not contribute 
In(p2) In(p'2) when k / p’, k' / p. The idea behind this is very simple: if momentum k / p’ 
enters the lower vertex in 17j, the fermion propagators are knocked “out of alignment” with 
p, and we lose the In(p”). This argument also establishes, e.g., that 17i does not contribute 
In(p?) In(p’?) in any momentum configuration. 


A.B. Carter, CH. Llewellyn Smith / Perturbative QCD 421 


Fig. 21. 


= Oe ae 
ry 


Fig. 22. 
p’ 
k ||P! 
A B 
Fig. 23. 
> —. 
| ( k II! | 
A B 
Fig. 24. 


tively: 
2 2 
24 Q Q 
wrt =—1¢h in( £5) in( $5) | 


The only other p’? logs in W,, are provided by the renormalization factor Z§ (p'2/u?): 


p? 1 p?\ | 11 4 p”? 
zi(Ex) =1+5A Cf in(25 + gd? Cp [Cr — (3 Co — 37 F)] in (2) : 


It is now a simple matter to verify that all the p’* singularities cancel in the final 
result. There remains a residuum of this cancellation in the form of a factor Z§(Q?). 
Thus the effect of radiation from the p’ line is to provide a factor Z§(Q7), just as 
corrections on the gluon lines produce Z3(f) (see CC6). This remarkable (though not 
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unexpected) result is represented by the blob in fig. 25. 

Final collection and comment: We have now shown to O(g*) that the leading 
logs can be assembled into ladder structures as in fig. 26. Let us remind ourselves 
what the notation stands for and obtain the final result. 

(i) The shaded circles represent u/traviolet renormalization factors only. Most of 
them combine conveniently to provide a running coupling constant at each rung, 
since (see fig. 27) 


[Z3())? Z3()_2O 
~TFF(AI2.g2* (3.12) 
[Z7()] 3 
(ii) The shaded half-moons on the side-rails correspond to exponentials like 
(3.11). Taken together, they provide an overall factor V, 


QO? 112-2742 
dk 2*(k¢) 

V(Q?, p?) = lac —. 

(Q p*) exp ah k? g? ’ 


(3.13) 


to W,. 

(iii) The remaining prescription is to evaluate fig. 28. The f, t’ integrals obviously 
lead to an exponential just like V, but containing 7,, instead of v. The corrections to 
the electromagnetic vertices satisfy a Ward identity which shows that each is equal 
to [Z£(Q?/u)]~!. Thus one e.m. vertex correction cancels the blob on the top line 
and we are left with [Z§(Q?/u2)]—! multiplied by Z§(p?/y?) from wave-function 
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renormalization (indicated symbolically by blobs on the incoming lines). The final 
answer is 


Q* 412 (22/42 
dkt (2°(kf) 
Wr = eX EE =| 
n Pp n—-1 B k? 2 
Ss 2 ax z°(x) NS P* ax 7 (x) 
_ N 
conf nts fF (LO) one, f S(E2)), on 
Pl n if x \ g n 1 x\ g ) 


which leads directly to eq. (1.3). 

Eq. (3.14) is the correct answer, and fig. 28 (with up to infinitely many rungs) 
the correct prescription, to all orders, although at this point we have only checked 
it to fourth order. In subsequent sections we justify our prescription to all orders. 
We proceed explicitly, although we could have worked inductively if we had been 
prepared to assume infrared and KLN finiteness: in each order only a few new ele- 
ments are introduced and finiteness is sufficient to determine them uniquely (it is 
easy to see that our results for parallel logs due to fig. 17 could have been derived in 
this way). 

As far as we know, ours is the first explicit calculation of W,, to O(g*). It can be 
used to deduce #?(t) to O(g*) or, equivalently, the function B(g) to O(g3). This is 
interesting since heretofore 6 has been deduced from “unphysical” quantities. Writ- 
ing 

Eee ae Se ee 
Br (t)=gr +A 6m ers ; 


various diagrams contribute to A as follows: 
A(11a + 15a + 15c) =—3Cg, 


A(12a, UV) =Cg—2Cp, 
A(12b, UV) =—3Cg, 
A(12c) =26;, 

A(8a + 9a) =5Tr, 


A(total) =—CgtiTp. 
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Note that: 

(i) The contributions proportional to Cy are the same as in the abelian case and 
cancel due to the abelian Ward identity Z, = Zp. 

(ii) The screening contribution proportional to Ty, due to qq production, is also 
essentially the same as in the abelian case (where it gives the total contribution). 

(iii) The remaining contributions are essentially non-abelian in character, being 
proportional to Cg. The negative (antiscreening) part is due entirely to diagrams 
involving the trighuon coupling: the ultraviolet contribution of fig. 12b and the paral- 
lel logs of 12b and 12k combined in 15A (we regard the role of wave-function renor- 
malization, 15c, as being merely to change the scale of 11a + 15c from k? to x”). We 
can, if we wish, interpret the decrease of #7(r) caused by these diagrams as being 
due to a “spreading out”’ of the fermion colour charge seen by the bifurcated gluon. 


4. Parallel logs and Ward identities 


In this section we derive the parallel log factor (eq. (3.13)) and obtain the corre- 
sponding factor for the case of external gluons, which is needed in the mixing sec- 
tion. We consider explicitly the factors associated with incoming particles; the fac- 
tors for outgoing particles are obviously the same. 

Consider fig. 29, where here and below the symbol! 1 means that the amplitude 
is one particle irreducible in the channel composed of the particles it brackets (the 
incident quark channel in this case). Otherwise the amplitude is complete, i.e., it 
contains all corrections to the other external lines. Apart from the incident quark 
and virtual photon, all the other external particles are physical — on-shell and trans- 
versely polarized in the case of gluons. 

We proceed up the incident gluon line until we encounter the first “‘parallel”’ 
gluon which we exhibit in fig. 30a. Note that this must be associated with a k / p 
log: if k were parallel to an outgoing line p + & would not be parallel to p and we 
would lose a log. Dropping terms which do not produce singularities as p? > 0, we 
can replace the factor (f — &) Y, from the quark propagator and vertex by 
2(y/(1 — v)) Ky, as in eq. (2.9). The factor k, allows us to use Ward identities. Since 
these relate complete amplitudes it is convenient to replace fig. 30a by 30b. 

The relevant identity is shown in fig. 31, where the dashed line is a Faddeev- 
Popov ghost field (¢%), a and i are colour indices, J/g and J}, are sources and @ is the 
gauge-fixing parameter (or gauge-fixing function a(k*) in our case). When the ghost 


Fig. 29, 
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Fig. 30. 
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Fig. 32. 


passes through the blob it can interact in all possible ways (note that the sources of 
the other particles could also serve as ghost “sinks” if they were not physical). The 
interactions with the sources are J§d" Ai, and WJi on the left (suppressing spinor 
indices) and J§¢* and g¢* w;4r 5Aj)J{ on the right. This is a special case of the Slavnov- 
Taylor-Ward identities derived by" t Hooft and Veltman [14]. It expresses directly 
the invariance of (0|T(¢* AA ... Aw ... W)|O) under Becchi-Rouet-Stora (BRS) trans- 
formations [15]. 

Using fig. 31 we can rewrite part of 30b in the way shown in fig. 32 (there are 
other terms, generated by both 32a and b, in which the ghost couples to the inci- 
dent quark source but they cancel). The first term (fig. 32c) on the right does not 
contribute to singularities as p? > 0 and can be dropped in LLA. We denote its con- 
tribution by 


d*k 1 
nr Seeruee 
Sk*k* p-Kk-&£ 

On the Jeft we have 


1 
aM, > 


MK 
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where M,, represents the total (one particle irreducible) contribution of the upper 
parts of 32a and b. However, we have used (p — £) Y, to generate k,, so we actually 
need 


1 em ee - (_* — a (p — K) + (32d term) 

ke" (p—kP 4 kPR2(p-kP (p-k-#) 
The 32c term can therefore be dropped in leading log approximation since (p — k) ap 
in the parallel configuration and # kills the potential log from the p line if it is vir- 
tual (or gives pu (p) > 0 if it is real). 

Thus only 32d survives. In this diagram the ghost quark interaction involves only 
the single scale (p — k)* ~ k? =k?. Therefore, it only contains UV logs (with argu- 
ment k?/u2) and we may redraw the bottom half of 32d as shown in fig. 33. Finally, 
the parallel log due to 30a may be represented by fig. 34 where the diagram on the 
left has the meaning discussed earlier (fig. 7) and the blob on the right represents the 
function of k?/p? due to @(k) X all UV logs. We now collect and analyse the latter 
factors which are: 

(i) a(k?), introduced by the Ward identity; this is just the contribution of the 
correction to the bare gluon propagator (see eq. (1.4)) in fig. 30a; 

(ii) The corrections to: (a) the quark gluon vertex in fig. 30a, (b) the quark propa- 
gator in fig. 33, (c) the ghost propagator in fig. 33, (d) the ghost quark Ward identity 
vertex in fig. 33. 

The factors (iia) and (iib) above combine with \/a(k?) to turn one factor of g 
into #(k?) in the lowest-order prescription for the parallel log factor. The behaviour 
of the remaining factors is determined by their anomalous dimensions, which are 
obtained from the factors which renormalize them. These are: Z3!/? for V/a(k?), 
Z3! for factor (iic) and Z 1Z3//Z, for factor (iid) (this follows from the fact that 
Zo is renormalized by Z 3/2/Z, and, as shown by Brandt [15], Z, AP" W; is a finite 
operator if ¢ and y are renormalized fields). The product of these factors is 

Z1Z3 

Z3Z, 
which is just one by the Slavnov-Taylor identity [17]. The net effect is, therefore, 
that the other factor of g in the lowest-order parallel log factor should also be 


replaced by 2(K?). 
We can now proceed up the incident line in fig. 34 treating subsequent parallel 
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Fig. 36. 


gluons in the same way. The limits of the k? integrals are determined by the neigh- 
bouring invariants, which set the scale of deviations from the simple (k?)~! factors 
in the integrands. Thus, for example, corresponding to fig. 35 we obtain 


Kf 


t n2 2 «Kk. 2 2 
_9 7,09, dk aes dk ae dk 
p2 t p2 t p t 


in leading log approximation. Finally summing all the parallel factors we obtain the 
exponential factor V of eq. (3.13). 

It is worth remarking that this technique can be used to sum leading logs which 
are contributed by crossed ladders in the case that we choose to perform covariant 
spin sums for produced gluons (which means that we must also include ghost pro- 
duction). This leads to effective ladder diagrams (as discussed by Gribov and Lipa- 
tov [3] in the much simpler abelian case). 

We also need parallel log factors associated with external gluons as illustrated in 
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fig. 36 (we already encountered such logs in (3.6)). The argument presented for the 
quark case goes through with only minor modifications. The result is that with 

g? > 27(k?) eq. (3.6) describes the factor in fig. 36b and summation produces the 
factor represented in fig. 36c: 


1 2 /32¢L2 
2-y 1ty dk¢ ( &*(kz) 
0 t 


5. The complete non-singlet result to leading order in Z* 


We have now established to all orders all the elements of the prescription for cal- 
culating leading logs presented at the end of sect. 3. The related prescription for 
calculating “fragmentation functions” in LLA can be obtained by the same argu- 
ments; it is illustrated in fig. 37 where we have opened up the blob on the top rung 
of a ladder to obtain a rainbow. 

This leading log prescription cannot be used directly to obtain realistic results 
since the leading log expansion of Z7(f) ~ 1/In t is only adequate for large t and the 
invariants in our diagrams have lower limits of order p?. The cure involves the intro- 
duction of another scale M, large enough that we are prepared to neglect #7(M?) 
compared to 1, which is used to separate parts of the diagram carrying momenta 
smaller and larger than M (this procedure is the diagrammatic realisation of eq. 
(1.2)). 

Our prescription is useful wherever the invariants are large and implies that only 
ladders can sustain large momenta (f) to leading order in 8*(t). Therefore to leading 
order in g*(Q) every diagram begins as a ladder, as shown in fig. 38 line 1 (where 
we omit UV blobs and parallel log correction factors for simplicity). We can split the 
t-integral into the regions |¢| >M? and |t!< M7, as indicated in line 2. For It|<_M? 
we can say nothing since ¢ runs down to O(p?). Accordingly we draw a “‘black box”, 
representing an unknown function (x, ft) (V for non-singlet) which gives the proba- 
bility of finding a quark with “mass” equal to ¢ (\t| <.M@*) and momentum fraction 
x. For |t| >M? large momentum enters the next vertex which must therefore have 
a ladder structure to leading order in Z*(M7). Proceeding to split up the r’ and sub- 
sequent integrals in the same way we arrive at the last line of fig. 38. 

We now write the parallel log factor (see (3.13)) at the bottom rung as V(ty, ft) = 
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2 2 
Iti <m*e |t] > M 
= 2 
> 
Ith<m2 , ‘ se 
[ee] 
= ty 
2 
to It] >M 
n=-0 ty 
[t| < m2 
Fig. 38. 


V(tn, M*) V(M?, t) and associate the second factor with the black box. Integrating 
over ¢ and taking moments, our ignorance of the hadron vertex is isolated in 


1 M2 
Ny(M*) = { x1 dx f < V(M?, t) N(x, t). 
0 p2 


For M? = Q? = |t\|max only the first term in line 2 survives and we see that Np,(M?) = 
W,,(M°). For Q? > M? this is modified by the ladder sum for scattering on a quark 


of “mass” M: 


— Or ae ml dt 
[Z§(Q?)]! [= (A Crrn) id, 73 B7(t1) oi = em) 
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Q* ay =2772 NS ,,2 
a _|& (Q*) | 27 n— 1/8*8o 
xexo [cru f . roo)-{E 2 , 


where the provenance of each factor should by now be clear. The factorization of 
the ladder sum into a function of Q? and a function of M* ensures that the result is 
independent of the choice of the large scale M, which is necessary for consistency. 
The final result is, of course, eq. (1.3). 


6. Singlet structure functions and mixing 


The singlet structure function, like the non-singlet, is computed by exposing lad- 
ders wherever large momentum transfer flows, as in fig. 39. There are now two pos- 
sibilities for the black box: the hadron vertex can deliver up a gluon or a (singlet) 
quark. The relative probabilities for these two processes to occur is absorbed into 
(incalculable) factors G,(M7) and S,,(M), respectively, analogous to N,(M7) in 
sect. 5. Our task is to evaluate the Q? behaviour of the ladder sums. (The blobs at 
the e.m. vertex and on the top rung combine to give a factor [Z§(Q?)]—}, as before.) 

The only new feature in the singlet case is the possibility of having gluon side- 
rails in the ladder, leading to new rung elements shown in fig. 40. Fig. 40a is already 
familiar to us: the q > qg skeleton gives a factor ACrry, the blobs combine (using 
(3.12)) to give g*(t), and the half-moons signify the virtual radiation factor (3.13): 


dt dx _ 
An~ CrIn~ 2 (t) exp| A Cr vf S700] ; 


et] <M? > 
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Fig. 41. 


The (infrared-finite) skeleton factors associated with 40b,c are easily calculated 
(see the appendix). The blobs in these elements do not quite combine to give #7(t), 
the difference being ratios of Z factors: 


Z t 
Bn~ 99 SP BO eh cer f EHC], 


t 


d dx 
Cr ~ 189, OTE xed Cov [FHC]. 


The (infrared-divergent) skeleton factor for 40d is obtained from (3.3) by inserting 
y”, giving a factor p, defined in the appendix: 


dt_, * dx _, 
Dn~® Copn-- #0) exp} Cev f —#*(x)] . 


At first sight it seems that the factors of Z3 and Z§ in 40b,c might prevent us 
obtaining the expected prescription. However, these factors combine in a simple 
way when we connect 40b and 40c with pure gluon ladders constructed from 40d. 
The portion of a ladder shown in fig. 41a contributes: 


7 seg dt Zi) dt Z3(1) 
fig. 41a |—78S Bt) 12 sO [- WH BO t ol 


x Pa Cr Pn) f a Ut ee f =S Pm) 
mt t 1 t m 
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r’ 


x exp| A Cov f =| ~ exp|A Ca» f” = A| (6.1) 


ty t 


ae" t ax >, . at 
= |-799 z(t’) | exp Ros + 2yp) f es c9| [189 a) “|, 
t 


where we have expressed (p,, + v), Z3, and Z§ in terms of y,4 and yp, using the for- 
mulas in the appendix. Likewise, 40c and 40b are connected by pure quark ladders 
constructed from 40a, giving for fig. 41b: 


» dat’ 
fig. 41b ~ 8 PO) S| 


x exp 88 +298) f See] [-e9e0"]. (6.2) 
t 


It follows from (6.1), (6.2) that, if we intend to perform the ladder sums, we can 
associate with the vertices the effective matrix factor 
( B 


mp) ~ Un 2e FPO, (63) 


where I’,, is the anomalous dimension matrix defined in the appendix and J the unit 
matrix. The elements of matrix (1, — 2yr/)” give the factors associated with all 
m-rung ladders (with incoming and outgoing quark or gluon lines) which can be 
constructed out of figs. 40a—d. The Q? dependence of the ladder sum of fig. 39 is 
then given by the eigenvalues d* — 2y¢ of the matrix T — 2yp/ (d* are the eigen- 
values of I’). The extra dependence involving 2y¢ is removed by the blobs on the 
top rung, giving the final result that W,, is a linear combination of two terms with 
Q? dependences [g?(Q?)]24" 8760, 


7. Concluding remarks 


Our methods provide explicit verification of factorization, etc., in LLA and 
allow us to derive the ladder structure, which is the basis of the parton model (and 
the Kogut-Susskind [18], Altarelli-Parisi [12] interpretation of the RGE), using ele- 
mentary techniques — much of the physics being apparent in O(g?). Beyond the 
leading log approximation the virtues of our approach decrease rapidly. It can be 
used to demonstrate factorization to all orders but the terminology and machinery 
needed become increasingly cumbersome and this is not compensated by the ability 
to calculate explicitly or by new insights. All-orders proofs seem simplest in the 
axial gauge (although here also the identification of hard scattering kernels, anoma- 
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lous dimension, etc., with diagrammatic structures is complex beyond leading order 
and does not offer much insight). 

Information about jet structure, multi-jet cross sections, etc., in deep inelastic 
scattering is easily obtained using our methods. For example, it is clear that the 
dominant contribution consists of two jets and that there is a subdominant O(a,) 
three jet structure, with a cross section given by the prescription of Ellis et al. [19] 
and Floratos [20] (this is shown explicitly in ref. [13] where earlier arguments [21] 
are extended to the non-abelian case). 

Our techniques have many other applications. For example, they can be used to 
rederive and interpret predictions for deep inelastic scattering on a photon [22], 
which were originally derived using the OPE. Furthermore it should be clear that 
they apply to other hard processes as well as deep inelastic scattering (e.g., to the 
production of heavy y-pairs). 

As an example we consider the production of heavy quarks (Q) in pp collisions, 
where the proper use of QCD perturbation theory is especially easy to understand 
using the diagrammatic approach. Many phenomenologists have assumed that fig. 42 
gives the dominant contributions, where the open circles represent the momentum- 
dependent quark or gluon distributions and the subprocesses is evaluated using g at 
some author dependent scale. This is correct in LLA in the (impossible!) limit s > °% 
with Ma /s fixed (typical choices of g only differ in the next to leading terms); 
in this limit 


2 
~fryin = M+ 45 (1 —|/1 -) +o, §>4M2>@~, 
so the large f[§] in 42a[42b,c] arranges the rest of the diagram into ladders, as in 
deep inelastic scattering [heavy u-pair production]. As s > % with Me fixed, how- 
ever, it is clear that QCD perturbation theory will become useless since —f pi, > 0. 
The picture is only useful, therefore, for Ms large and § “‘not too far from thresh- 
old” (in practice (§) may remain of order 4M, up to quite large s due to the steep- 
ness of the quark and gluon distributions). 
Quantitatively we obtain precise predictions for the moments of the cross section 
corresponding to fig. 42. These momenta factorize and the predictions have the 


: Q s s FS 
q a 
B Cc 


Fig. 42. 


2 


(=) 
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(symbolic) parton model interpretation 


oo 


1 1 
f o(s, MR) Se= fi x""q(x, MQ) dx fx" q'(x', MB) dx 
0 0 


1 


ea names 
Xf 6G Ms. 
1 


where z = s/4M, z= §/4MB. The integrals only converge for n 22 (orn > 1 in the 
non-singlet case), so these precise predictions are indeed only sensitive to s near 
threshold, as anticipated and interpreted above. 

This is just one example of the utility of the diagrammatic interpretation of hard 
scattering processes. 


We are grateful to John Taylor for numerous useful discussions. One of us 
(A.B.C.) thanks the Rhodes trust for the award of a Rhodes Scholarship and the 
other thanks the S.R.C. for support. 


Appendix 


In this appendix we compile the notations and conventions used throughout this 
paper. Rules for Feynman diagrams are those found in ref. [23], with the exception 
that our spinors are normalized so that uu = 2m. We use the standard convention 
that QO? = —q? = q? — q2 >0 in deep inelastic scattering. The effective expansion 
parameter is 


§ 


2 
Sr2- 


~ 
Ill 


A. 1. Figures 


In figures a wavy line ((“VV\V) is a gluon, a dashed line (— — —) a ghost, a solid 
line ( ) a quark, the inductance symbol (eeeee) a virtual photon, and a dotted 
line (----: ) indicates the effective contribution (summed up with Ward identities) 
of virtual parallel radiation. 


A. 2. Lie group factors 
Cy = quadratic Casimir eigenvalue for fundamental (quark) representation, 


= — 1) for SU(N) , 
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Cg = quadratic Casimir eigenvalue for adjoint (gluon) representation 
= N for SU(N) , 


Tr= 3 (number of flavours). 
A.3. The renormalization group and anomalous dimensions 


The normalization of our anomalous dimensions and 6-function is fixed relative 
to renormalization group equation (written here symbolically; G is some function 
with definite renormalization group transformation properties): 


7) 3 
[2 +82-1]e-0, (A.1) 


This normalization may differ (by factors of two) from those used by other authors. 
Here 


B(g) = — Bog? + O(g') , 
with 
1 
Bo =o [$Co —- <Tr] 


The running coupling constant 2(Q7) satisfies 


og(Q?)_ _ 
adin(Q2w) FP? 
and is approximately 
g 


OG cdi sain ee oe waa 
BO T+ Fogg? ins?) 


The anomalous dimension matrix of twist-two operators in the operator product 
expansion, appropriate to the nth moment of W(x) is 


ee 2S, ) 


rT ae od 
i 8 8S 


—yNs, = ACESn ’ 


fr=-3]1 a ar 


oe + 
(n+ 1Xn+2) j-2 J 


2+3n+4 
_vyQG =9T. Fespeeerercy : 
aioe Flint 1X(n + 2)(n + 3) 
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2 
_.6Q = n*+3n+4 | 
rad Mee +1)\(n+2)J’ 


+1 
i 4 4 " 
GG .! | > 1, 4 
—ys™ =erxCnrl—- + + = 4 ~(— 2r[sT, 
TA GLB (nt 1) (nt2Xn+3) a I 2X[37 Fl. 


Our Z factors are defined relative to renormalized full propagators and one par- 
ticle irreducible vertices normalized to their bare value at euclidean points charac- 
terized by mass yw: 


p°_@—— = 26 (p2/u2) (————] 
p?---@ -- = J; (p2/p2) (- ------ ) 
p2 p2 
p2 
p* p? 
IPI 7 
ZF (p2/p2) 
p2 
28 W2 ‘ ’ 
P ie 2 . 1 ‘a UU” 
ipl. © SS : (A.2) 
9 ¢ Z, (p*/*) 
p 


neglecting induced terms in the vertices with different tensor properties (which are 
negligible asymptotically) and setting m = 0. 

In our stagnant Feynman gauge, these functions satisfy renormalization group 
equations like (A.1) with the appropriate field anomalous dimensions expressed in 
the stagnant Feynman gauge (which are the same as in ordinary Feynman gauge in 
O(g2)) and without the usual “gauge drift” term 63/da. Thus, e.g.,Z3(P?/) satis- 
fies 

) 


—d a 
Lames "Pag ~ 270] Za) =O; = 
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_ 1 5 4 
Ya =—4rAl-6Co + 37 FI 
and solution 


47 4/8789 


; 


Z3(P2/u 2) = &(P") a 


ZE satisfies (A.3) with yf = —4rCp, etc. 
The Z’s satisfy the identity 
[Z3(p?)1? _ [23(e7)]? Z3(p?) _ 8°(e”) 
Z1(p?)]?— (ZF (e))? B 
and the Slavnov-Taylor identities 
Z3(p?) _Z5(p?) _ Z3(p) 
Zi(p2) ZF(p2) Z,(p?) 


A.4. Infrared-divergent radiation factors 
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(A.4) 


(A.5) 


(A.6) 


Associated with real gluon radiation from an incoming quark line is the factor 


pre: tn = 


ro is appropriate to an outgoing quark line. 


Associated with virtual gluon radiation (summed up with Ward identities) from 


quark lines (incoming or outgoing) is 


E ue fo 


Associated with real gluons from an incoming gluon line is 


a 


(Po for outgoing gluons) and for virtual radiation 


Ps 


Real and virtual combine to form the infrared-finite factors 


y y 
+y(1—y) 


ty 
= 


tn tvs — Yn Si +2rF], 
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At this point we also list the finite factors (used in sect. 6) associated with 
(i) quark emission from a gluon: 


pPntv= (vES, + 27,4] - 


2n? + 6n+8 


1 
| J 9" 490 1 Ee yen + ayn 3) 


—1 
= ,QG . 


(ii) quark emission from a quark: 


2i-—y)  n*+3nt+4 


1 
: | dyy” yt 
poy 4 n(n + 1\(n + 2) 
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